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1. Introduction 

Classical, non-dissipative solutions - called "lumps" or "defects" - have been studied exten- 
sively in the history of physics with various motivations behind their investigation. Initially, 
they were studied in the context of Scott Russell's waves in hydrodynamics , then as "soli- 
tons" by mathematical physicists and more recently by particle physicists, as possible new 
"particles" in the spectrum of non-linear field theories. They have been investigated by 
condensed matter physicists looking into superconductivity and by astrophysicists study- 
ing the formation of galaxies. Even without considering the fiuids and condensed matter 
research, a vast amount of literature ^'"^'4'^'" has accummulated on different aspects of such 
solutions and the interest continues unabated. 

Lumps have been observed in fiuids, plasmas and a number of condensed matter sys- 
tems but, in particle physics, not a single solution of this type has been detected. The 
magnetic monopole remains elusive and, in fact, the archetypal GUT monopole is perceived 
to be enough of a cosmological problem that it has to be "infiated" away. Cosmic strings, 
though astrophysically promising, have yet to be detected and are tightly constrained by 
the millisecond pulsar and other observations. Heavy domain walls are believed to be 
cosmological disasters and a particle physics model is considered inadmissible if it predicts 
them. The question still looms large if there are any classical, non-dissipative solutions in 
particle physics. 



On the other hand, the detection of a defect in a system can give us valuable information 
about the system. Since the defect is non-perturbative, it gives us information about the 
non-perturbative structure of the theory. The existence of a topological defect would tell 
us something about the topology of the theory and with it, other features that would be 
impossible to glean by perturbative scattering experiments. A famous example is that due 
to Dirac : the very existence of a monopole would tell us that electric charge is quantized*. 
The presence of lumps in a system can lead to novel phenomenon: GUT monopoles can 
catalyze proton decay ', and, strings and textures can lead to galaxy formation and to 
the generation of primordial magnetic fields iu,ii,i2,i3 L^]2lpg can also give rise to exotic 
quantum phenomena such as fermionic zero modes and quantum hair on black holes . 
The benefits that would be reaped if lumps exist in particle physics seem to far outweigh 
the doubts one may have about their existence. It is hardly surprising, then, that so much 
effort has gone in the past several decades to uncover the mysteries of the lump. 

In the following notes, we shall merely touch upon certain basic aspects of this immense 
subject and hope that this can be a starting point for the reader to follow up on the 
references that have been provided. The choice of topics included here are essentially 
aspects of the subject that I have been personally involved in together with some basic 
underlying topics that are the foundations on which the subject has grown. 
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It should also be said that this reason is no longer as compelling as it used to be in view 
of the fact that electric charge is automatically quantized within the framework of Grand 
Unified theories. 



2. Topological defects: field theory 

In this section we will study topological defects as classical solutions in certain field 
theories . 

The general criterion for the existence of ad (spacetime) dimensional topological defect 
in a field theory which exhibits spontaneous symmetry breaking from a group G to a 
subgroup H is: 

n:,_d{G/H) + 1 (2.1) 

where TXn{G/H) is the n homotopy group of the coset space G/H. The cases (i = 0, 1, 2, 3 
correspond to the texture, monopole, string and domain wall. The condition (2.1) leads 
to topological defects but the theory might also contain non-topological, semilocal and 
embedded defects whose existence cannot be detected by using (2.1). 

It is often more convenient to think in terms of the vacuum manifold, Ey, described 
by the values of the scalar field (f) that minimize the potential: 

r dV d^V ^ , , 

If the surface Ey has incontractible surfaces of 3 — d dimensions, the theory will contain 
a d spacetime dimensional topological defect. 

A physical justification for the above criterion can be given. Suppose that the con- 
figuration of (/) on a 3 — (i dimensional surface at spatial infinity (5") is denoted by (^^o. 
If we assume that there is vanishing energy at infinity, then (/)oo lies on Yjy. Therefore 
(^oo describes a mapping from S to Ey. Next let us imagine the case when the image 
of this mapping is one of the incontractible ?> — d dimensional surfaces in Ey and denote 
this image by /. But 5" is contractible - assuming that space itself has trivial topology* 
- and so we can continuously shrink this surface to a point. When 5" shrinks to a point, 
it must continue to be mapped to a non-trivial surface in Ey since / is assumed to be 
incontractible. But this would mean that (^ would be multi- valued at the point to which 5" 
has been contracted and this is not acceptable. The only way out of this contradiction is 
that the field (/> must leave Ey at some point in space. However, this means that (^ cannot 
remain at the minimum of the potential everywhere and there must necessarily be at least 



If there are black holes or other gravitational peculiarities present, these arguments would 
need modification. 



one point where there is non-zero potential energy. The location of this potential energy 
is the location of the topological defect and the energy distribution at this location defines 
the energy distribution of the defect. Note that this argument shows that the asymptotic 
configuration of the field (j) is sufficient to determine the existence of the defect. 

The case of the texture is somewhat different since the asymptotic field configuration 
is not sufficient to determine its presence and at almost all times the field never leaves 
the vacuum manifold. There is energy in the configuration because the symmetries are 
global and hence the variations in the field carry gradient energy. The topology can be 
understood by considering the configuration in spacetime since then we can consider three 
spheres in the vacuum manifold that are incontractible. If the model one is considering 
only contains gauged symmetries, there are no textures since all the gradients in the scalar 
field can be compensated by gauge fields. Yet one may still have textures in the entire 

universe and these result in different sectors in the gauge theories and lead to degenerate 

4 
vacuua . 

To make things more concrete, we now list the simplest models that give rise to walls, 
strings, monopoles and textures. 

The domain wall solution arises in models in which a discrete symmetry is spontaeously 
broken. For example, consider the model: 
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{d^<pf - A(</>2 - rj''f\ (2.3) 

where, (^ is a real scalar field. The symmetry breaking in this model is Z2 ^ 1 and hence 
txq{G/H) 7^ 1. In terms of the vacuum manifold, it is given hj (f) = ±r] and consists of two 
disconnected minima. The consequence is a domain wall in the model which interpolates 
between the two minima. For a static domain wall in the yz plane, the field solution is: 

(j) = ?7tanh(vA rjx) (2-4) 

and the energy per unit area of a wall in the y^— plane is: 

Eyj = ^VXv^ . (2.5) 

The most familiar example of a model with gauge strings is the Abelian-Higgs model: 

1„ „,„. A/,0 v^\^ 



Ss = j d^x \{d^ + ieA^)ct>\' - \f^.F^''' - ^(</>2 - "Ly 



(2.6) 



where, (p is a, complex scalar field and F^i/ = d^Ajy — djyAfj^. Here the symmetry breaking 
is U{1) -^ 1 and niiG/H) ^ 1. The vacuum manifold is given hj (j) = ^e*" where a is 
any phase. Therefore the vacuum manifold is a circle (parametrized by a) and contains 
incontractible 1 dimensional curves. The corresponding unit winding string solution in 
cylindrical coordinates (r, 6, z), and along the 2— axis, is of the form: 

</> = ^/(r)e^^ (2.7) 

A^ = -^d^9 (2.8) 

er 

where, the functions /(r) and v{r) vanish at the origin and go to 1 as r ^ oo. There is no 
known closed form for / and v and they have to be found numerically. 

When 8A = e , that is, when the scalar and vector masses are equal, the energy of the 
vortex can be found analytically using Bogomolnyi's method and the result is: 

Es = TTTj^ . (2.9) 

For other values of the parameters, the energy has to be evaluated numerically . 
A simple example of the magnetic monopole ^^'^^ occurs in the model: 

1„„ „.„. A/^ ?72\21 



Jd^x \{d, + zeeM-)^? _ ^_F-^F-^^- -^{^-y)' 



(2.10) 



where, (/> is a triplet of fields, a = 1,2,3, F/^j^ are the non-Abelian field strengths and 
{€°')ij = eaij is the usual epsilon symbols. The symmetry breaking is 0(3) -^ 0{2) and 
this gives magnetic monopoles since Tr2{G/H) ^ 1. Here the minimum of the potential is 
a two sphere and hence contains incontractible two dimensional surfaces. The monopole 
configuration can be written down in spherical coordinates: 
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(2.11) 



At^ = ^a^Jr'^,f^[^^) (2.12) 

where, 

f = {sin9 cos(j)j sin9 sincj), cos9). (2-13) 

The functions /(r) and v{r) and the mass of the monopole need to be found numerically. 
The only exception is in the Prasad-Sommerfeld limit, when the solutions and the mass 



of the monopole are known in closed form. This is the hmit A -^ and with e and ry held 
fixed. Then the solution is: 

/(r) = Crcoth{Cr) - 1 (2.14) 

Cr 

<r) = — ^ (2.15) 

sinh[Cr) 

where, C = rje. The energy in this limit is: 

E = '^. (2.16) 

e 

Finally, the simplest texture ' occurs in a model with symmetry breaking 0(4) -^ 
0(3): 

St= I d^x [(a„$)2 - A($2 _ ^2)21 (2.17) 



where, $ is a column vector of 4 real fields. The vacuum manifold in this case is a three 
sphere and admits incontractible three dimensional surfaces. The texture solution is time 
dependent. The usual approach is to assume spherical symmetry and only consider the 
cr— model limit when |$| = r]. Then, in spherical coordinates , 

$ = r]{cosXj sinxr) (2-18) 

where, f is the unit radial vector as given in (2.13). The function x(t, r) can also be found 
analytically: 

X(t, r) = 2tan^^{-u) , u<0 (2.19a) 

X(t,r) =2ta?i"^(+w) +7r , l>w>0 (2.196) 

X{t,r) =2tan'^{+l/u)+n , u>l (2.19c) 

where, u = r/t. 

This completes our synopsis of topological defect solutions. The description is far 
from complete but suffices for applications to cosmology. The reader interested in the 
classification of defects, defects in particle physics models, exotic defects etc. is referred 
to the excellent review by Preskill and to the book by Vilenkin and Shellard . 



3. Embedded defects: field theory 

Even if the general criterion for the existence of a topological defect (eq. (2.1)) is not 
satisfied, the model can still permit the existence of topological defect like solutions ^"'^' . 
These solutions are essentially the topological defects of a smaller theory which are em- 
bedded in the bigger theory under consideration. The conditions under which such an 
embedding can be successfully carried out are not very stringent and so we can expect 
embedded defects to exist in almost any model. 

The existence of a solution in a model does not automatically mean that it is stable and 
it is in this crucial way that embedded defects differ from their topological counterparts. A 
single topological defect is stabilized by topology and is separated from the zero topological 
defect sector by an infinite energy barrier. Embedded defects, however, need not be stable. 
In fact, mostly they are unstable, sometimes they are metastable and the only known 
examples of stable embedded defects are semilocal strings . 

Instead of giving the general arguments for constructing embedded defect solutions, 
we shall only present some illustrative examples. We shall first construct an embedded 
domain wall solution as this is the simplest example of an embedded defect and then we will 
focus on the electroweak model and construct the embedded electroweak string solutions. 
A limiting case of electroweak strings will give us the semilocal string. 

Walls - The most trivial embedded solution is a domain wall embedded in a global 
G = U{1) model. We express the Higgs field in terms of two real scalar fields (p"', a = 1, 2. 
A Lagrangian that is invariant under the global U{1) rotation and describes static field 
configurations is, 

L = a,(/)«aV°-A ((/)''(/)''- rj'^] , (3.1) 

with i labeling the spatial coordinates. 

The first step in constructing the embedded domain wall solution is to identify a Z2 
subgroup of the full symmetry group. Let us consider the Z2 subgroup defined by the 
transformation: {(pi,(p2) ~^ {~4'1t4'2)- Any non-zero vacuum expectation value of (^1 will 
break this Z2 subgroup completely and so the embedded symmetry breaking is Z2 ^ 1. 
This symmetry breaking has topological domain walls: 

(pl = ?7tanh(vA r]x) (3.2) 



and so this configuration for (pi togetlier witfi (^2 = is our candidate embedded domain 
wall solution. 

Once we have identified a candidate embedded defect solution, we should check if 
it extremizes the energy functional. The general conditions for this to be true can be 
written 26,27 ^^^ require setting up a formalism. The idea however is simple: to check that 
the configuration is a solution, we perturb the configuration and verify that the variation in 
the energy vanishes to first order in the perturbations. Since we know that the domain wall 
is a solution to the theory when (p2 is zero, and the directions of (pi and (p2 are orthogonal, 
there is no need to perturb (pi - only perturbations in (p2 might be dangerous. Now we see 
that (p2 appears quadratically in the energy functional and so the variation in the energy 
functional vanishes to linear order and the configuration is a solution. 

This argument for checking when embedded configurations are solutions can be ex- 
tended to arbitrary models and defects without much difficulty. 

Electroweak strings '^ - 

Consider the Weinberg-Salam ' model of the electroweak interactions. The sym- 
metry breaking is: SU{2)]^ x U{1)y -^ U{1) and the bosonic sector of the Lagrangian 
is: 

L = -\G^,aG^-- - \fb^,F^^'' + II^A'/'P - A(</>t</> - 7^2/2)2 (3.3) 

where, (/> is a complex doublet. The definitions of the field strengths are: 

Gl. = d^W^ - d.W^ + ge^^'W^^W^ (3.4) 

F^^„ = d^B^ - d^B^ (3.5) 

and the covariant derivative is: 

D, = d,-i^-T-W--ii^B, (3.6) 

where, g and g' are coupling constants and t"' are the Pauli spin matrices. 
The electroweak energy functional follows from the Lagrangian (3.6): 

\FB^JFB^J + (I^,</>)"^(I^,</>) + A(</>t</> - r,^ /2f (3.7) 



E = I dz I (fix 






where, i,j = 1,2,3 and we have restricted ourselves to the case when there is no time 
dependence and the time components of all gauge fields vanish. In addition, since we 



will only be interested in string solutions, we will only consider configurations that do not 
depend on the 2;— direction. Then the integration over z can be ignored and we can think 
in terms of the energy per unit length of the string. 

The first step is to choose a t/(l) subgroup of the full symmetry group. We choose 
this to be the t/(l) subgroup generated by 

r^ = -cos^e^r^ + sin^Oyjl = diag{- cos 2^^;, 1) . (3.8) 

(Note that T is the generator corresponding to the Z— boson (see (3.11) below)). 
Now the candidate embedded string solution may be written down: 

<l>emb = ^/.or(r)e*^'Vo- (3-9) 



where, we take. 



(3.10) 



Here, (r, 9) are polar coordinates. In addition, we want that the covariant derivatives 
vanish at infinity and so we take 

Zi^ = cos9wW^ - sin9wBn = [A^j^or (3-11) 

where, [A^]i;or is defined by eq. (2.8). All the other fields in the model are taken to vanish. 

One can check that the static configuration (3.9), (3.11) extremizes the energy func- 
tional and hence is a solution . 

Different choices of the ("embedded") subgroup lead to other string solutions. The 
choice that we now consider is the subgroup sitting entirely in the SU{2) factor of the 
electroweak model and generated by: T^ = sina r + cosa r , where a is some constant. 
Then the corresponding embedded string solution is: 

(Pemb = /«or(r)e'^"Vo , sinaWl + cosaWf = iA,)yor , (3.12) 

and all other orthogonal combinations of gauge fields vanish. 

The one parameter family of string solutions in (3.12) is called the VF(a) string since 
the fiux in the string is purely in the SU{2) sector. Furthermore, by a global gauge 
transformation, any single string solution in the family - that is, a string with any value 



of a - may be transformed into the string configuration with a = 0. Exphcitly, this gauge 
transformation is: 



i) = exp 



.(1 + r 



3^ 

(3.13) 



together with a corresponding transformation of the gauge fields. This does not, however, 
mean that if there are many different string solutions of different a present, they can be 
gauge transformed to another multi-string configuration with all strings having the same 
value of a. The simplest way to see that a is a non-trivial parameter is to consider a 
loop of W string such that a runs from to 27r as we go around the loop. The winding 
of a around the loop is a discrete number and cannot be altered by any non-singular 
gauge transformation. Hence, a loop with varying a is not gauge equivalent to one with a 
constant value of a. 

One may also see that strings of different a are distinct by comparing the directions 

of the field strengths in group space in each of the strings. Of course, the field strengths 

are only gauge covariant and not gauge invariant and so one must first parallel transport 

the gauge field of one string to the location of the other string and then take the scalar 

product of the field strengths. This leads us to the following quantity: 

/ r /'•^2 -1 \ 

A = Trfr"Fg(f2;a)P exp(-z / dl ■ W^t'') T''F[^{xi;a)) (3.14) 

The quantity A is a gauge invariant measure of differences in a between strings. 

Semilocal strings - 

Consider the case ^r = in the electroweak model. Now the symmetry group is 
SU{2)gi X U{1)y (where gl stands for "global") and it breaks down to a U{l)gi group. 
The vacuum manifold is still given by the minima of the potential and is a three sphere. 
However, the model continues to have the electroweak strings as solutions. But the Z-string 
is the only gauge string since there is only one gauge field in the model. 

With the knowledge of embedded defects, the existence of the semilocal string solution 
is not a surprise, but what is surprising is its stability. The simplest case where one 
can explicitly check the stability is when A = e . In this particular case, the method of 
Bogomolnyi can be used and it is at once obvious that the string minimizes the energy . 
A more careful analysis ' reveals that the string is only neutrally stable in this case, 
unstable for larger A and stable for smaller A. Subsequent numerical studies ' have 
confirmed this result. 



One way to understand the stability of the semilocal string is by inspecting the pro- 
cesses by which the string solution can destabilize. These processes necessarily require the 
presence of gradients of the Higgs field for which there can be no compensating gauge fields. 
Hence, unwinding requires a growth of the gradient energy but accomplishes a decrease of 
the potential energy. When the coupling constant A is large, the potential energy is the 
more important piece in the energy functional and the string prefers to unwind. If A is 
small, however, the gradient energy required to unwind the string is prohibitive and so the 
string is stable. 

Topological aspects of semilocal defects have been investigated in Ref. 38, 39 and such 
defects have been constructed in a wide range of theories in Ref. 40. Some cosmo logical 
aspects of semilocal strings have been investigated in Ref. 41. 

Next, by considering the case of small but non-zero values of g, it is clear that even 
the Z— string will be metastable for some values of parameters. A plot of the region of 
parameter space where the string is stable may be found in Ref. 42. The stability of 
electroweak strings (and other embedded defects) can be considerably enhanced if there 
are bound states present on the string . 



4. Cosmological formation of defects 

So far we have only discussed topological defects as classical solutions in certain field 
theories. What relevance can such solutions have for cosmology? Here we will argue that 
these defects would form during phase transitions in the early universe and, in most 
cases, would have survived until the present epoch. 

Consider a cosmological phase transition occurring at some temperature T. During 
the phase transition the Higgs field (j) acquires a vacuum expectation value so that \(j)\ = t]. 
If the vacuum manifold is non-trivial, this still does not fix the location of the vacuum 
expectation value on the vacuum manifold. For example, if the model is the one with 
Z2 symmetry (eqn. (2.3)), we can have (p = +1] or (/> = —77 at any given point in space. 
Furthermore, the value that is acquired at one point is uncorrelated with the value acquired 
at some other point provided the points are separated by a distance greater than the 
correlation distance ^ at the phase transition. Causality necessarily implies that ^ < t 
where t is the epoch of the phase transition. Hence, after the phase transition, the Higgs 
field lies on different points of the vacuum manifold at different spatial points. Then, it 
will happen, just by chance, that the configuration of the field on some asymptotic surface 
will be topologically non-trivial. (For example, in the Abelian-Higgs model (2.6), there 
will be closed circuits in space on which the phase of the Higgs field varies from to 2n as 
each of the circuits is traversed.) When this happens, a defect will necessarily be present 
somewhere inside the surface. 

This existence proof of the formation of defects during cosmological phase transitions 
relies on two facts: (1) the presence of a defect can be determined solely by looking at the 
asymptotic field configuration and, (2) the vacuum expectation value is uncorrelated on 
distances larger than t. This mechanism for the formation of defects is called the "Kibble 
mechanism" . 

Realistically, the phase transition is a thermal process and one should study the for- 
mation of defects using statistical field theory. Less ambitiously, one should estimate ^ by 
finding the thermal correlation function during the phase transition and this has been at- 
tempted in simple models. A rigorous estimate of, say, the density of defects after a phase 
transition is not available. However, for cosmological purposes this is not very relevant 
either. The point is that even if some defects - as few as one per horizon - are produced, 
the cosmological effects can be very dramatic. 



Another mechanism by which defects can be produced is by quantum mechanical nu- 
cleation of the defect during an inflationary phase of the universe . A quick way to 
understand this phenomenon is by thinking of inflation as a strong force that pulls apart 
any object. On the other hand, the mutual attraction of monopoles to anti-monopoles, the 
tension in cosmic strings and domain walls, tends to collapse these objects. In terms of a 
potential, there are two regions - one when there is no defect and the other where there is 
a defect-anti-defect pair (or loop of string, or shell of domain wall) that are being pulled 
apart due to inflation. These two regions are minima of an "effective potential" and they 
are separated by a potential barrier. But then there can be quantum mechanical tunneling 
from one region to another and if there are no defects to start with, inflation can literally 
pull a pair out from the vacuum. As this is a quantum mechanical effect, the number den- 
sity of defects produced in this way is small. However, in the post-inflationary universe, 
our horizon is only a small patch of the universe and it may well be that our horizon is 
located in the region which does contain a few defects. If these defects are domain walls 
or strings, they may still be relevant for cosmology within our universe . 

What do the defects look like on formation? The following summary is based on nu- 
merical simulations to study the formation of topological defects during phase transitions. 

Walls: There is one infinite wall and very few smaller walls whose size distribution is 
exponentially suppressed '. The fraction of wall energy that resides in the infinite wall 
is about 87% . 

Strings: There are a few infinite strings with close to 80% of the total string length . 
The remaining string is in loops with a scale invariant distribution: 

dn{R) = c-^ 

where, dn{R) is the number density of loops having size between R and R + dR and 
c ~ 6 is a numerical coefficient determined by simulations . The loops will collapse, 
radiate energy and disappear but the infinite strings will survive as they are protected by 
topology. 

Monopoles: The number density of monopoles is ~ ^~ where ^ is the correlation 
distance after the phase transition . For a second order phase transition ^ ~ T and 
for a first order phase transition it is given by the typical size of the bubbles when they 
collide. The radius of bubbles at collision could be anywhere between T and the Hubble 
distance ~ t. 



Embedded defects: The Kibble mechanism does not directly apply to the formation of 
embedded defects since the point (1) above is not met and, in this case, a full study of the 
phase transition appears necessary. However, some guesses can be made on the basis of 
what is known about the formation of topological defects. For example, if the embedded 
strings are metastable they can end on monopoles. And if the probability for a string 
to break by the formation of monopole and antimonopole pair is small, one would still 
expect to form infinite strings. For a larger probability, however, the length distribution 
of embedded string segments is expected to be exponential ^^'^^ - that is, the number of 
long strings is exponentially suppressed. At this time, however, no one has a quantitative 
understanding of the formation of non-topological defects. 

Nucleated defects: For defects that nucleate in de Sitter space, it is possible to find the 
wave-function for quantum fiuctuations of the defect and hence the distribution of various 
shapes and sizes of nucleated defects '. Cosmogical consequences of such defects have 
also been investigated in Ref. 53. 



5. Cosmological constraints on domain walls and monopoles 

Domain walls and magnetic monopoles are strongly constrained by cosmology. Consider 
domain walls first . 

As discussed in the previous section, an infinite domain wall will be formed during a 
domain wall forming phase transition. The infinite domain wall will move under its own 
tension and try to straighten out. Immediately after the phase transition, the motion of the 
wall is damped by friction but as the plasma gets diluted by Hubble expansion, the drag 
decreases and eventually the motion of the wall is effectively undamped by friction. Hubble 
expansion is still important. The single domain wall in the universe cannot disappear since 
it is protected by topology and would be present in the universe if it were ever produced. 
(Inflation could, however, push the domain wall outside our horizon in which case it would 
be irrelevant for our observable universe.) Assuming that the domain wall straightens out 
completely, its area within our horizon is ~ t and its mass is at where a is the energy 
per unit area of the wall. Therefore the energy density in the wall is ~ a/t and, for the 
walls not to dominate the universe today, we require that a/t < p, where, p is the matter 
density at time t. The energy per unit area of a domain wall is usually (for example, in 
the model (2.3)) given by a = vAry . Assuming that the coupling constant A is of order 1, 
and taking p ~ 3/(327rG't ) the constraint on domain walls gives rj < lOMeV. 

The scale ry is the scale at which the phase transition occurs and hence the only per- 
missible walls are the ones that can form relatively late in the history of the universe. This 
excludes domain walls that form at the Grand Unification scale or even at the electroweak 
scale. 

Next consider magnetic monopoles. There are several interesting bounds on the number 
density of monopoles - each involving different physics. Here we will only consider the 
cosmological bound which in itself is quite severe ''. 

The cosmological bound is that magnetic monopoles should not overdose the universe. 
This means that the energy density in monopoles pm today should be less than the critical 
density of the present universe. The smallest possible number density of monopoles at 
formation (rif) is one per horizon at that epoch: 

n; ~ -3 (5.1) 

7 



After formation they quickly become non-relativistic and from then on their energy density 
decays hke that of matter and redshifts as a(t)~ where a{t) is the scale factor of the 
universe. The energy density of the radiation dominated universe, however, redshifts 
faster, in proportion to a(t)^ . Therefore the ratio of monopole energy density to critical 
density (O^) at time t is: 

„,„(*) . ^^ (B.2) 

Pc{tf) a{tf) 

where, m is the mass of the monopole. With pc = 3/{32nGt ), a{t) oc t ' and m ^ Tj 

where, Tj- is the temperature at the phase transition, we find that Qm becomes one at a 

time t* given by: 

10^^GeV^8 



U ~ 10-*^. (5.3) 

V m / 

For GUT scale monopoles {m ~ 10 GeV)^ this epoch occurs at 10^ s - well before the 
matter era. (This justifies our use of a{t) oc t ' .) Therefore, if monopoles were formed 
at the GUT epoch, their energy density would completely overwhelm our universe and 
would overdose it. This is not observed and so some way has to be found to resolve this 
"monopole overabundance problem"*. 

The monopole overabundance problem is head-on in conflict with the philosophy of 
conventional GUTs which is that the electroweak and strong forces are unified in a Grand 
Unified symmetry group - a simply connected group - at an energy scale of about 10 GeV 
and that there is no new force that comes into play between the electroweak (10 GeV) and 
the Grand Unification scale. Coupled with standard cosmology, this philosophy implies 
that the Grand Unified symmetry group broke down to a subgroup with a hypercharge 
U{1) factor at the Grand Unification scale. But then condition (2.1) is satisfied for the case 
of monopoles. This means that within conventional GUTs, the monopole overabundance 
problem must be confronted. 

One obvious solution is to relax conventional GUTs and allow for the possibility of a 
U{1) factor in the Grand Unified group. Then monopoles will never form and there will 
be no overabundance problem. 

If one is unwilling to relax one's philosophy of Grand Unification, there is still a cos- 
mological solution and another particle-physics solution to the monopole overabundance 



The problem is even more severe when one considers the various other constraints on the 
present monopole fiux - such as the Parker bound or the bound coming from neutron 

stars . 



problem. The cosmological solution is to have an inflationary phase during or after the 
formation of monopoles . The inflationary phase simply dilutes the monopoles until their 
number density gets so small that they cause no problem. The particle physics solution 
is due to Langacker and Pi who consider the formation of monopoles at the GUT stage 
and then another stage when the electromagnetic symmetry gets broken. At this stage 
the magnetic flux of the monopoles gets conflned and the monopoles get connected by 
strings. The strings bring the monopoles and antimonopoles together, they annihilate, 
and subsequently the electromagnetic symmetry is restored. 



6. Cosmic strings: general properties 

Here we will summarize some properties of cosmic strings. 

At formation, the string network consists of closed loops and infinite strings. It is 
estimated that ~ 80% of the energy in the string network resides in infinite strings . 

Once formed, the strings move under their own tension and try to straighten out. This 
motion is damped due to the frictional force of the ambient matter ^"'^^ and is also slowed 
due to the Hubble expansion. The frictional force is more important than the Hubble 
expansion as the matter density is high. With time, however, the matter density gets 
redshifted and the Hubble expansion dominates the frictional force. The time at which 
the Hubble expansion drag become comparable to that due to friction is t^ ~ (Gn) tpi- 
After this time, the frictional force can be ignored . 

The motion of a string in vacuum with energy density ^ is well described by the 
Nambu-Goto action: 



S ^-ix j dTda\J-gy^) , (6.1) 

where, g^ ' is the determinant of the world-sheet metric defined by 

^6^ = g^,u^aX^'{T, a)d^x%T, a) , (6.2) 

a,h = T, a; g^u is the metric of the background spacetime and x^{t, a) are the coordinates 
of the string world-sheet. 

The Nambu-Goto action is valid as long as the radius of curvature of the string is 
much larger than the thickness of the string ''. Also, the action does not include 
the interactions of strings when they intersect. When there is a collision of a string with 
another string, it leads to the phenomenon of intercommuting ' ' . In this event, the 
strings intersect, exchange partners and then again move as given by (6.1). 

Another factor that plays an important role in the evolution of the string network is 
the gradual loss of energy from strings into forms of radiation. An oscillating string is a 
time dependent solution to a set of field equations and so one would expect that the motion 
would lead to the radiation of quanta of any fields that couple with the fields that form the 
string. However, it has been shown that the radiation in these quanta is negligible from the 
strings that are expected to be of cosmological interest . This is solely due to the fact that 
the curvature and oscillation frequencies of such strings is very small compared to particle 



physics scales and hence the only radiation that can possibly be emitted is into massless 
particles. A more detailed study then shows that the dominant loss from oscillating gauge 
strings is to gravitational radiation '. A loop of size L emits gravitational radiation 
and loses all its energy in a time: 

where the coefficient F is numerically found to be of order 100 for certain family of string 
loops . (In the case of global strings, the energy loss is dominated by the emission of 
Goldstone bosons '.) 

At times later than t,,,, the evolution of the network of strings is governed by string 
tension, Hubble expansion, intercommuting and gravitational radiation. These four factors 
make the evolution complicated enough so that no one has a clear picture of what the 
network looks like at any instant. Progress in this problem has relied on the results of 
numerical simulations of the string evolution and, recently, an analytical attack is also 
underway. 

Even though the string network is not fully understood, a few features seem to be 
emerging '^J-Jj '4,75,76,77 pj^-g^ Qf g^^^ ^j^g network at any instant much later than t^ seems 
to consist of a few infinite strings (that is, strings that traverse the whole horizon without 
closing up on themselves) and a large number of tiny loops. The size distribution of the 
loops is not known but the favoured guess is that the size is given by the gravitational 
radiation cut-off distance: / ~ TG^t. (Loops smaller than this evaporate in less than 
a Hubble time and would probably not be significant for cosmo logical purposes.) The 
distribution of loops in space is not known either but, since the loops are produced at very 
high velocities (f ~ 1) one would expect them to be distributed roughly homogenously 
even if they are initially produced in a localized region of space (for which there is some 
visual evidence). The long strings are not smooth but have a lot of irregularities. The 
scale of the irregularities is guessed to be the same as the size of the loops ~ TGfxt at any 
time t. These irregularities are called "kinks" or "wiggles" in the literature and the long 
strings are said to be "wiggly". 

Cosmic strings can also have the ability to carry persistent electric currents . Such 
superconducting cosmic strings can have very dramatic cosmological signatures. Other 
notable varieties of strings include global, non-abelian and Alice strings. We shall not 
discuss these varieties of strings but the reader can find a description in Ref. 6. 



7. Cosmic strings: gravitational properties 

Non-superconducting topological defects interact with their environment primarily via 
gravitational forces. Here we will consider the metric of gauge strings and textures and 
describe some of the known properties . 

We first consider the metric of a source with energy-momentum tensor °^'°l 

T;; = 6{x)S{y)diag{i^,0,0,T). (7.1) 

With T = n = hq this is the effective energy-momentum tensor of an unperturbed string 

with string tension ^q as seen from distances much larger than the thickness of the string. 
When ^2,83 

I^T = IJ.I (7.2) 

this also describes the energy-momentum tensor of a wiggly string as seen by an observer 
who cannot resolve the wiggles on the string. 

The gravitational field of the string can be found by solving the linearized Einstein 

^00 = %3 = 4G'(^ - T)ln{r/rQ), 



equations with T^ from (7.1). This gives 



hn = h22=4G{^ + T)ln{r/rQ), (7.3) 

where, h^i, = g^u — rj^i, is the metric perturbation, r = {x + y ) ' and tq is a constant 
of integration. 

For an unperturbed string, T = fx = hq and we get 

^00 = ^33 = 0, hii = h22 = 8G^oln{r/ro). (7.4) 

A coordinate transformation brings this metric to a locally flat form, 

ds'^ = dt^ - dz'^ - dr'^ - (1 - 8G//o)r^#^ (7.5) 

It describes a conical space, which is just a Euclidean space with a wedge of angular size 
Aq = SttG/Uo removed and the two faces of the wedge identified. A particle at rest with 
respect to a straight string experiences no gravitational force, but if the string moves with 
velocity Vs, then nearby matter gets a boost 

Ui = ATvGfXQVs^s (7.6) 



in the direction of the surface swept out by the string. Here, 7^ = (1 — Vg)^ ' . This 
effect is responsible for the formation of wakes and for a discontinuous change of the 
microwave background temperature across a moving string '. Assuming that the string 
is perpendicular to the line of sight, the magnitude of the latter effect is 

— = 87rG^QVs-fs ■ (7.7) 

The conical metric also results in the formation of double images of background objects. 
In the cosmological context, this would lead to the gravitational lensing of background 
quasars and galaxies '. 

Returning now to the wiggly string metric (7.3), we first consider the effect of the 
wiggles on light propagation . Assuming for simplicity that the direction of propagation 
is perpendicular to the string, we can write the relevant components of the metric in the 
form 

ds^ = (1 + hoo) [dt'^ - {dx^ + dy^)] (7.8) 

where we should identify the half-lines y = ±4nGnx, x > 0. The conformal factor (l + Zigo) 
does not affect light propagation and can be dropped. Then the resulting metric describes 
Minkowski space with a deficit angle SttG/U, and we conclude that background temperature 
discontinuities produced by wiggly strings are given by the same equation (7.7) with fXQ 
replaced by fx. In contrast to the smooth string, however, the wiggly string also produces 
a change in the photon temperature when the photon propagates parallel to the string but 
perpendicular to the velocity of the string. 

Next, we study the formation of a wake behind a moving wiggly string. First look at 
the problem in the rest frame of the string where the particles are flowing past the string 
with a velocity Vg in the x-direction. The linearized geodesic equations in the metric (8) 
can be written as 

2x = -{l-x^-y^)d^hoo, (7.9) 

2y = -{l-x^~y^)dyhoo, (7.10) 

where over-dots denote derivatives with respect to t. We need only work to first order 

in Gfx, in which case (7.10) can be integrated over the unperturbed trajectory x = Vst, 

y = yQ. Then we can transform to the frame in which the string has a velocity Vs- The 

result for the velocity impulse in the y-direction after the string has passed by is : 

2nG{iJ-T) ^ ^ ,^^^, 

^^ = ^Jl L _ 4nGiJVs-fs (7.11) 

Vsis 



The second term is the usual velocity impulse due to the conical deficit angle. But, for 
small velocities, it is the first term that dominates the defiection of particles. The origin 
of this term can be easily understood. From eqn. (7.3), the gravitational force on a 
non-relativistic particle of mass to is F = 2mG{iJi — T)/r. A particle with an impact 
parameter r is exposed to this force for a time At ~ r/vs and the resulting velocity is 
u, ~ {F/m)At ~ G{p - T)/vs. 

The metric of a loop of cosmic string can be obtained quite easily in the weak field 
approximation '. The result is that, at distances much larger than the size of the 
loop, the metric is Schwarzschild with mass parameter M equal to the mass of the loop. 
Singular points ( "cusps" ) on the string world-sheet may have novel gravitational features 
but a proper treatment of these features requires going beyond the weak gravitational field 
approximation. (In addition, these singular features do not seem to be very relevant 
for cosmology since now it is believed that loops are not very important and that the 
occurrence of cusps on loops is not as generic as it first seemed to be°°.) 



8. Structure formation by wiggly cosmic strings 

The formation and evolution of long-string wakes and their possible role in structure for- 
mation in the scenario where strings are not wiggly have been discussed in Refs.84, 89, 90, 
91, 92, 93 and many other papers. The scenario where the strings are wiggly has a shorter 
history and some of the relevant papers can be found in Refs. 13, 94, 95, 96, 97, 98. Here 
we will follow Ref. 13. 

When a coUisionless fluid flows in the wiggly string metric, the gravitational field 
focusses the fluid particles inwards so that streamlines flowing on either side of the string 
converge behind the string. Therefore a wake forms behind the string which has twice the 
density of the ambient fluid. If the fluid is not coUisionless, the fluid flow into the wake 
will be accompanied by the formation of shocks and turbulence. These features of the 
wake are likely to be important for the formation of structure on galactic scales and for 
the generation of magnetic fields. However, the details of the wake are not important for 
the formation of large-scale structure which is what we shall now describe. 

Consider a wake formed behing a moving string segment of length ~ C{ti) at time t^. 
The distance travelled by the string in one Hubble time is ~ Vgti, and thus the initial length 
and width of the wake are Ij ^ tj, Wi ^ Vgti. We shall first assume that the universe is 
dominated by cold dissipationless matter. In this case the two opposite streams of matter 
in the wake overlap, and the mass density is enhanced by a factor of 2 within a wedge 
with an opening angle 2ui/vs, where Ui is from eqn. (7.11). The average thickness of this 
wedge is di ~ -u^tj. The initial surface density of the wake is 

a, « 2p{t,)d, ^ -^^ (8.1) 

3 Vsti 

where p{t) = {6nGt ) is the average density of the universe, and its total mass is 
Mi ^ CTihwi ^ (n — T)ti. Note that M^ is independent of the string velocity Vg- If the 
string moves faster, the wake is wider, but the surface density is decreased proportionately. 
We note also that the velocity perturbation (7.11) is produced at distances up to ~ Wi 
from the plane of the wake. For r > Wi, the gravitational field of the string is like that of 
a stationary rod and ui ~ G{fx — T)ti/r. 

As the universe expands, the length and width of the wake grow like the scale factor, 
a{t) ~ t ' , while the total mass of the wake grows by gravitational instability like M oc 



a{t). As a result, the wake thickness (defined as the turnaround distance) and surface 
density evolve like d oc a (t), a oc a^ (t). At the present time {t = Iq) 

ao^'-^iff. (8.2) 

Vsto \tij 

Cold-dark- matter wakes can also be formed during the radiation era (t^ < teq), but in this 
case the gravitational instability sets in only at t ~ tpq. It can be shown that the surface 
density of the resulting wakes is proportional to (ti/teq) ' . Together with eqn. (8.2) this 
implies that the most prominent wakes having the largest surface density are the ones 
formed at t ~ teq- 

The fraction of the total mass of the universe accreted onto wakes which were formed 
at time ~ ti can be estimated as (for ti > teq) 

f - f||y - 87rG(/. - T)z, (8.3) 

where Zj is the redshift at t^. The total mass of dark matter in all wakes is dominated by 
the wakes formed at t ~ teq, 

hot ~ 'i^G^lQZeq ~ 0.4/i2^6- (8.4) 

Here, h is the Hubble constant in units of 100/cm/s.Mpc, the universe is assumed to have 
critical density, O = 1, (Ug = G/Uq/IO^ , and in the last step we have used the values of n 
and T from the simulations. 

The evolution of the initial velocity perturbation (7.11) can be found from the equation 
of motion for dark matter particles, 

d 

u-\ — u = g (8.5) 

a 

where, g = 2'KGa{t) is the gravitational acceleration due to the wake. This gives u{t) oc 
t ' . A careful analysis shows that the present velocity perturbation due to a single string 
impulse is 

2 1/2 /o «N 

UQ ^ -UizJ . (8.6) 

If the filamentary wakes created by strings were sheet-like, they would have a thickness 
~ Uitizf. For wakes produced around teq-, this thickness is somewhat larger than the width 



of the wake ~ VstjZi and hence we should not treat the filamentary wake as having planar 
geometry. Instead the wake should be treated as having linear geometry and the accreting 
structure will be cylindrical in shape with the possibility of some planar sub-structure. 
The diameter of the cylindrical structure is characterized by the geometric mean of the 
previously calculated width and thickness and is ~ (uiVg) ' z^ tj while the length is tiZj. 

To get a qualitative feel for the appearance of the wakes, we adopt the picture developed 
in Ref. 72 for the evolution of the string network. The basic idea is that the long strings are 
moving slowly at speeds ~ 0.2 for about one Hubble time. Then there is an intercommuting 
somewhere in the network which triggers an instability and speeds up the string to a much 
higher velocity ~ 0.6. In this way, during every Hubble time period, a string moves slowly 
for most of the time but the slow motion is followed by a rapid motion that helps maintain 
the scaling solution in which the distance between strings stays a fixed fraction of the 
horizon size. In addition, string simulations show '-^''^ that the coherence length of strings, 
beyond which the directions along the string are uncorrelated, is ^(t) ~ t. The inter-string 
separation L{t) is of the same order of magnitude. In the matter era L{t) ~ 0.7t. The 
rms string velocity on the scale of the smallest wiggles is (< v >) ' ~ 0.6, but the 
coherent velocity obtained by averaging over a scale ^ is Vg ~ 0.15. The average mass per 
unit length and string tension are (in the matter era) fx ~ 1-4/Uo, T ^ 0.7fXQ. With these 
values, the first term in eqn. (7.11) is about ten times larger than the second. 

If a string segment moves coherently for more than one Hubble time, the resulting 
wake will have a variable surface density, with denser parts being the ones formed at 
earlier times. The straightening of long strings on the scale ^ ~ t occurs mainly due to 
string intersections. Long strings occasionally self- intersect producing a horizon-size loop 
which then rapidly collapses into miriads of tiny stable loops. If two different strings 
intercommute, the highly curved regions near the points of intercommuting develop a high 
velocity, Vg ^ 1, and also shed off a large number of tiny loops as they move. The wakes 
due to rapidly moving strings have the form of sheets with dimensions tiZi x t^Zi x Uitizf 
while the wakes due to slow strings have a filamentary appearence. As we explained (see 
below (8.1)), the masses of both types of wakes are comparable, but the surface density in 
the filamentary wakes is much higher, and we expect filamentary features to be prominent 
in the large- scale galaxy distribution. In addition to wakes due to long strings, there will 
also be comet-like wakes produced by rapidly moving small loops . The characteristic 
scale of the large-scale structure in this scenario is teqZeq ~ lOh Mpc. With h = 0.5 it is 



comparable to the scale suggested by observations^^ (~ 25h Mpc). 

The wiggliness of the string network implies that the wakes will not be uniform but 
will have sub-structure on the scale of the wiggles. This scale is expected to be larger than 
the damping scale due to gravitational radiation from the string network which is TGiit, 
at the time of formation of the wake, where F ~ 10 is a numerical factor coming from the 
rate of gravitational radiation"^. For the wakes produced at tgg, the sub-structure is on 
a comoving scale larger than ~ I^qK^ kpc. We expect that the wakes will fragment into 
smaller objects due to this sub-structure. 



The large-scale velocities predicted at the present time can be found from eqn. (8.6). 
For sheet-like wakes from rapidly moving strings, it gives uq ~ 300/Ug h km/ s where we 
assumed that t^ ~ teq and Vs'^s ~ 1- These velocity perturbations extend over regions of 
size {lQh~ Mpc) and may account for the observed large-scale streaming velocities . 
Reasonable values of uq are obtained, e.g., for h ~ 0.5, /Ug ~ 4. We note that in some 
regions of space the motion of matter can be affected by two or more different strings. 
The streaming velocity in such regions will typically be enhanced by a factor ^/n where 
n is the number of string impulses that the matter experiences . Regions larger than 
{IQh Mpc) will also get peculiar velocities due to string impulses but the velocity will 
scale as 1/L where L is the size of the region. (This is simply because a string gives a 
coherent impulse to a region of size L, when its correlation length becomes comparable 
to L. For large L, this happens later, giving less time for the velocity to grow.) The 
observational situation on the dependence of peculiar velocity on length scale is quite 
unclear and it remains to be seen if this predicted fall-off agrees with observation . 

If the dark matter is cold, wakes formed at all epochs prior to radiation-matter equality 
will survive and density fluctuations will be present on very small scales too. Albrecht and 
Stebbins ' argue that this small scale power is excessive and the sheet-like structures 
formed later would not be prominent. The situation to me does not seem as clear since 
one could imagine small scale structures themselves clustering into larger scale structures. 
So the large scale wake could be prominent simply because it rearranges the small scale 
structure into sheets and filaments. 

In a universe dominated by light neutrinos, wake perturbations are damped by neutrino 
free streaming on co-moving scales smaller than Aj/(t) ~ Vi,{t)t, where Vi,{t) ~ Veq{t(>q/t) ''^ 
is the rms velocity of neutrinos and Veq ~ 0.2. On larger scales the evolution of pertur- 



bations is similar to that in cold dark matter. For a cold-dark-matter wake formed at 
time ti, all matter initially within a distance u{tizi/ z will be accreted onto the wake by 
the redshift z. A neutrino wake will go nonlinear at the redshift Zr^i when the co-moving 
scale of \y{tj) becomes less than the distance to which the matter has been swept by the 
wake : Wjtj(a^//aj) ~ X,^{ti)a^i/ai, where the scale factor a{t) is related to the redshift 
hy 1 + z = a(to)/a(t). For filamentary wakes, this gives '^ 1 + 2^; ~ A.bfiQh indepen- 
dent of Zi- With h = 0.5 and fXQ = 4, we have 2„; ~ 3.5. For sheet-like wakes, we find 
1 + Zj^l ^ 2/jQh . Observations do indicate that 2 = 2 — 3 is the epoch of intensive galaxy 
and quasar formation . Thin wakes of small relativistic loops are strongly suppressed 
by the neutrino free streaming , and it appears that loops play a negligible role in this 
scenario. Eqn. (8.4) then implies that most of the matter in the universe remains unclus- 
tered at the present time . This may explain why dynamical measurements in clusters 
give values of O substantially smaller than 1. By contrast, in the cold dark matter scenario 
the loops accrete at least as much matter as the wakes, and the voids will be pierced by 
the long comet-like wakes formed behind relativistically moving loops. The characteristic 
scale of the large-scale structure in this scenario is teqZeq ~ lOh Mpc. With h = 0.5 it is 

comparable to the scale suggested by observations (~ 2^h Mpc). The surface density 

-1/3 
of the neutrino wakes produced subsequently decreases but the decrease is only oc t^ . 

This means that the structure on still larger scales can also be prominent. 

Baryonic wakes in a neutrino-dominated universe start collapsing after baryons de- 
couple from radiation, t > t^^ec- However, since baryons constitute only a small fraction 
of the total density, the growth of these wakes is strongly suppressed. Baryonic wakes 
could nonetheless be cosmologically significant if the energy output from the primordial 
stars formed in the wakes might trigger some kind of explosive amplification and lead to 
preferential galaxy formation along these wakes . They could also explain the existence 
of quasars at redshifts greater than 3. The scale of baryonic wakes, tf^ec-^dec ~ 50/i~ Mpc, 
is comparable to the largest-scale structure observed in the universe. 

A novel outcome of the cosmic string scenario is that it predicts the generation of 
primordial magnetic fields ^■^"°^^^', The mechanism by which this happens relies on the 
fact that the relativistic motion of strings after decoupling of matter and radiation induces 
vorticity in the baryonic fluid. The vorticity then leads to the generation of primordial 
magnetic fields. 

To put these arguments on a firmer footing it is necessary to establishe generation of 



vorticity and then to show that the vorticity wiU lead to magnetic fields. Fortunately, the 
second step had been investigated in the 70's and several mechanisms are known by which 
vorticity can lead to magnetic fields. So the main task that remains is to show that there 
will be vortical motion and to estimate the vorticity. 

An estimate of the Reynold's number for the flow of the baryonic fluid into the string 
wake shows that it is very high (~ 10 ) and hence it is natural to suspect that the flow will 
be turbulent. A large Reynold's number, however, is not sufficient to guarantee turbulence 
and one must demonstrate that the fiow is unstable. In the case of fluid flow into a cosmic 
string wake one can actually go further and explicitly describe the mechanism by which 
vorticity is generated. 

Consider the wake formed by a relativistically moving string at the recombination 
epoch. At this epoch, the sound speed is dropping steeply but the flow of the fluid into 
the wake is still given by (7.11). Then, for string tensions that are suitable for structure 
formation, the fluid flow is supersonic and the wake is bounded by strong shocks. On 
large-scales this shock is uniform but on small scales, the shock is non-uniform because the 
string is wiggly and the wiggles have highly variable velocities. This is the crucial feature 
- the wiggly string wakes are bounded by strong, non-uniform shocks. 

Once we have shown that the scenario has this feature, the presence of vorticity follows. 
From Euler's equation and steady flow, one flnds that the vorticity is related to the gradient 
of the entropy by the equation : 

vxuJ = -TVs (8.7) 

where, v is the fluid velocity, to the vorticity, T the temperature and s the speciflc entropy. 
In the preshock region, the flow is isentropic but, at the shock, the entropy suffers a 
discontinuous jump. Since the shock is non-uniform, the post-shock entropy is different at 
different points along the shock. This gives us gradients in the post-shock entropy and a 
non-zero vorticity. 

The presence of vorticity in the baryonic fluid flow means that the protons and electrons 
are in vortical motion. But we also have ambient photons and neutral particles which 
interact with the protons and electrons. The next crucial ingredient in the scenario is that 
the masses of the protons and electrons are different and therefore their interaction times 
with photons are also different. (The scattering cross-section of either particle with photons 



is inversely proportional to the square of the mass of the particle.) Then the stronger 
interaction of photons with electrons slows them down with respect to the protons and the 
resulting electric current due to the differential rotation of charges produces a magnetic 
field. Such scenarios - using the different interaction rates of protons and electrons with 
photons - were proposed by Harrison , Mishustin and Ruzmaikin and others in the 
70's. 

While the scenario is clear qualitatively, quantitative estimates are more difficult to 

obtain. The first hurdle is to understand the vorticity in the fiuid fiow. For this we 

would need to have a quantitative analysis of the turbulence in the fiuid fiow and, as far 

as I know, there is no theoretical recipe for analyzing turbulent fiow. But dimensional 

arguments allow us to estimate the average vorticity as follows. We expect that vorticity 

will be produced on the scale on which the wake is inhomogenous. Therefore the co-moving 

scale of this vorticity is ~ FGntiZi where Zj is the redshift at the epoch tj when the vorticity 

is generated at a time t^. The velocity of the fiuid is estimated as in (7.11) and so the 

vorticity is 

u 0.1 ,^ ^, 

W ~ -r ~ _ (8.8) 

Given the vorticity, the generated magnetic field is estimated from the results of Mishustin 
and Ruzmaikin 

B « 2!^4±il!^. (8.8) 

where, m is the electron mass, e the electron charge, z the red- shift at which the vortical 
motion starts, O the mean-to-critical density or the baryonic matter, H the Hubble con- 
stant, Te"f{0) the interaction time between electrons and photons at the present epoch and 
uj is the angular velocity of the eddy. The interaction time between electrons and photons 
at the present epoch is given by 

where, ay ~ 10^ cm is the Thompson scattering cross-section and ^7(0) is the present 
photon energy density. Inserting O ~ 0.03, if ~ 50 km/ s~Mpc, z = 10 , and lv = 10^ s 
gives, 

B ~ lO"^'^ G . (8.10) 



The magnetic field produced due to the vorticity at decouphng can be further amphfied 
by turbulence in the wake and by a galactic dynamo. Such a tiny seed field is all that is 
needed to generate the observed galactic magnetic field of 10^ G. 

Cosmic strings seem to be uniquely suited for generating magnetic fields via vorticity 
since they naturally have the two features that seem essential for this mechanism: coherence 
and strength. The first feature is that the vorticity should be on relatively large scales 
(~ lOkpc) for which it is essential that the source producing the vorticity should also have 
this coherence scale. If the source is not lineal , it is difficult to see how this large a scale 
is obtained. The second point is that the vorticity on these scales has to be relatively 
large, implying that the source itself has to be undergoing violent motion. Once again, 
this feature comes up naturally in the string picture but seems hard to get with other 
sources. 

Outlook: 

While the cosmic string scenario for structure formation has ingredients that seem to 
be promising, it is not detailed enough yet to be testable. This is because the problem 
is doubly difficult - first one has to understand the evolution of the string network and 
then the evolution of the matter that gravitates around the network. The evolution of the 
network itself has turned out to be a really sticky problem and an analytical understanding 
is just beginning to emerge. The fiow of matter around the string network promises to be 
an even more difficult problem since, as we have seen while discussing the generation of 
magnetic fields, the fiow will be non-linear and turbulent. Without going into the details 
of the fiow, one can only make some broad predictions about the large-scale structure as 
we have done above. To make predictions on galaxy scales, it is essential to understand 
the structure and fragmentation of the wake. Only then will we be able to say something 
about the galaxy-galaxy correlation function and other quantitative measures that can 
enable a comparison with observation. 

These problems seem so difficult that at times I am tempted to think that it may be 
simpler to observe cosmic strings or rule them out (as a mechanism for structure formation) 
on the grounds that they are not observed. The millisecond pulsar observations seem to be 
a foolproof way to go but for this we have to wait for another decade or so. The microwave 
background anisotropy measurement by COBE does not confirm or rule out strings. But 
small-scale measurements - when they become more reliable - could test the cosmic string 



scenario*. The direct observation of strings by their gravitational lensing property is also 
possible but is likely to be effort consuming. On the other hand, this effort seems very 
worthwhile considering all the exciting outcomes! 
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